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A proper version of the proto renormalization-group scheme is presented to derive ampU- 
tude equations in striped pattern formation with conserved and nonconserved order parame- 
ter. In the conserved case, the resuh preserves the conservation law as well as the rotational 
invariance of the physical system. This makes a great contrast with existing singular per- 
turbation theories such as the multiple-scales method, in which the perturbation expansion 
truncated at any finite order destroys those symmetries. 

Subject Index: 056, 054, 013 

§1. Introduction 

Spatial patterns are seen in a great variety of physico-chemical systems. They 
are typically periodic in space, at least locally. The most common are stripes (or 
often called rolls or lamellae) and hexagons. The unity of dynamical mechanisms of 
their birth and evolution allows us to consider some models of the nonlinear theory 
of such structures. The Swift-Hohenberg (SH) equationP is one of the simplest and 
most canonical paradigms, and has been intensively studied in the past.'^ 

In order to describe the appearance and growth of crystalline phases that occurs 
in many technologically important materials processing, a dynamical model called 
phase-field crystal (PFC) model is also introduced.l^ Since its introduction, the 
PFC model has been applied to a wide range of problems such as solidification, 
elastic and plastic deformations.'^ It is a continuum equation as the SH equation, 
and additionally there exists a global conservation law for the order parameter in 
this model. The presence of the conservation law is expected to give rise to new 
universal features of pattern formation which largely differ from those bestowed on 
the systems lacking the conserved quantity such as the SH model. ^' As a matter of 
fact, the PFC model is a conserving analogue of the SH model. 

Since it is impossible to obtain analytical solutions to those models in pattern 
formation problems that are nonlinear partial differential equations, it is quite natu- 
ral that a method to coarse-grain the model has been sought. Among others it is now 
relatively standard to use the so-called amplitude equation formalism. In the vicinity 
of the pattern forming instability, fluctuations of the order parameter field favor the 
formation of structures built on some finite wavevectors. The evolution of these un- 
stable modes is much slower than the evolution of the stable ones. Hence the faster 
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evolution of the stable modes may be adiabatically eliminated from the dynamics. 
As a result, the dynamics can be reduced to the slow dynamics of the amplitudes 
of these unstable modes. The amplitude equations that result have been a powerful 
tool in the study of pattern formation, selection and stability properties-l^ As such 
there exist many theoretical methods to derive them. Important examples include 
the method of multiple scales^' and renormalization-group (RG) methods.EJ® In 
particular, the proto-RG method which has been proposed in Ref. 8) is the most ab- 
stract version of the RG approach to free as much as possible the RG procedure from 
the necessity of obtaining explicit expressions of secular terms that are required of 
the other RG method.'^ It is therefore not surprising that the proto-RG scheme has 
been used to study the SH and PFC models .'^1' 'Ell For the present account, it is to be 
noted in the latter paper that the authors conjectured that one would need separate 
RG schemes to coarse-grain order parameter equations depending on whether or not 
the order parameter is conserved. Yet the issue remains controversial,'^ and in 
the present paper we have reformulated the proto-RG reduction of the SH and PFC 
models in a successful attempt to accommodate the dispute. 

In the present paper we will thus concern ourselves with a system with a real 
order-parameter field ■0 whose dynamics are described by the zero-noise Langevin 
equation 

dti^ = (id^r [eV' - V'' - (dl + k^M , (M) 

where p = for a nonconserved order parameter (NCOP), p = 2 for a conserved 
order parameter (COP). The dx is a two-dimensional gradient operator (V) with 
respect to the position vector x, and is the Laplacian (V^). Hence the choice 
p = corresponds to the SH equation, while the PFC equation is represented by 
p = 2. The e is the bifurcation parameter. It is well known that the competition 
between the surface energy contribution given by d^i^ term in the square brackets 
on the right-hand side (RHS) of Eq. (jl-ip and the curvature energy term (x d^i^ 
gives rise to spatially modulated structures with period ~ 27r//co for e > 0. In the 
present discussion we restrict ourselves to striped patterns for the sake of brevity, 
although the generalization to hexagonal patterns is straightforward. In doing so, 
we demonstrate in the following §2 that by a suitable modification of the original 
scheme,® proper amplitude equations can be derived on the same footing in the 
proto-RG formalism in both cases of NCOP and COP. As the essential logic of the 
formalism can be recognized at the lowest nontrivial order, in the present paper the 
explicit calculations will be exhibited to 0(e). We compare the analysis with the 
multiple-scales theory in §3, and concluding remarks are given in §4. 

§2. The Proto-RG Scheme 

Since we treat the term etp — i/j^ in Eq. p-ip as a perturbative one, we scale ip 
as ^/e^l^ and denote the new ip with the same symbol. The model (jl-ip then reads 

dti^ = {id^r KV' - V'') - (dl + fco') V] • (2-1) 

Throughout we take the average value of ip is zero. 



RG for Amplitude Equations for Cellular Patterns 



3 



We expand the solution to Eq. ()2-ip as 

^{t,x) = i^o + ei^i + e^iP2 + -- - ■ (2-2) 
The = 0, 1, 2, • • • , obey the equations 

L{^t,^^)^|Jo=0, (2-3) 
L{dt,d^)i;i= {id^f (Vo - V'o), (2-4) 
L{dt,d^)i;2= (iS^f {(l-3^g)Vi}, (2-5) 
etc., 

with 

m, d^) ^dt + {id^r {di + kif. (2-6) 

We have the zeroth-order solution from ()2-3p in the form 

^o(a;, B, A) = B{p) + (A(p)e*'^-" + c.c.) , |fc| = fco, (2-7) 

where B and A are, respectively, the real and complex function of an arbitrary space 
parameter p, which will be later used as the regularizing point. Here and hereafter 
we must set B = Q for the NCOP case, i.e., for p = 0. The dependence of the solution 
on B^A is explicitly denoted for clarity of the subsequent discussion. Thus 

x) = B+ (^yle^*'-^ + c.c.) + e^i{t, x, B, A) + e'^Mt, x,B,A) + -- - . (2-8) 

Notice that = e**''^ is the eigenfunction of the zero eigenvalue of L; L^q = 0. It 
is also the case with = constant when p = 2; henceforth we may arbitrarily set 
'Fq = 1. Hence and l^'o can be the source of the secular terms in > 1. 

Let ipj{r,R,B, A) be ipj{r,B^A) with variables r = {t, in the secular pref- 
actors of <Po and replaced by i? = {t, p}, discarding the constant terms. The 
regularization of ipj is made by introducing the renormalized B and A, denoted by 
Br and A/j, via 

A = ArZ = Ar{R) [1 + eZi{R) + e^Z2{R) + •••], 

B = BrW = Br{R) [1 + eWiiR) + e'^W2{R) + ■■■]■ (2-9) 

Then Eq. (j2-8p can be written as 



ip = Br% + {Ar<Po + c.c.) + e ^BrWiI'o + {ArZi<Pq + c.c.) + i:i{r, Br, Ar) 
BrW2% + {ArZ2^o + c.c.) + Mr, Br, Ar)+ 

BRWidBMr,BR,AR) + {ARZidAMr,BR,AR) + c.c.)] +■■■ . (2-10) 
The renormalization constants Zj and Wj are determined order by order via 

BrWi% + {ArZi^o + c.c.) + Mr, R, Cr) = 0, 
BrW2% + {ArZ2<Po + c.c.) + Mr, R, Cr) + 

BRWideMr, R, Cr) + (^ARZidAMr, R, Cr) + c.c.) = 0, 

(2-11) 
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where Cr = {Br, Ar}. It is then easy to see that Eqs. ()2-lip and ()2-10p are obtained 
by the e— expansion of the following equations, respectively: 

{B - Br)% + {{A - Ar)^o + c.c.) + e^i(r, R, C) + e^2{r, R,C) + ••• = (), (2-12) 



i/j = Br% + {Ar<Po + c.c.) + 



Mr,C)-iJiir,R,C) 



+ 



(2-13) 



In the above C = {B,A}. Notice that all the secular terms in ijjj's are removed in 
Eq. (|2T3]1 . 

Introduction of the variables R is equivalent to splitting the derivative dt to 
dt + dr, and dx to dx + dp. Since ip (and ipj) must be independent of the arbitrary 
regularization point R, we have Tip = and Tipj = 0, where T = L{dt + 9ri + 
dp) — L{dt,dx)- Applying the operator T to Eq. (|2-13p . we then find 

Hdt + dr, dx + dp) [Br% + {Ar^o + c.c.)] = eTMr, R, C) + e^TMr, R,C) + -- - . 

(2-14) 

Let V be the projection onto (pQ. Then by the definition of tpj we have 

VL{dt,dx)iPj{r,R,C) = 0. (2-15) 

Hence 

VL{dt + dr,dx + dp)AR^o = eVL{dt + dr,dx + dp)Mr, R, C) + 

e^VL{dt + dr,dx + dp)Mr,R,C) + -- - . (2-16) 

Similarly, using the projection operator onto ^q, Q, we have 

QL{dt + dr,dx + dp)BR% = eQL{dt + d^, dx + dp)Mr, R, C) + 

e^QL{dt + dr, dx + dp)Mr, R,C) + -- - . (2-17) 

Equations ()2-16p and ()2T7p are the basic equations of the proto-RG scheme. We 
remark that the argument of tpj is C, not Cr. Thus, in the proto-RG formulation, 
there is no such ordering ambiguity as suggested in Ref. ^ for renormalization and 
differential operations. Consequently, there should not be a different ordering (renor- 
malization before differentiation) as the authors proposed for the COP equation. 

We shall now carry out the 0(e) computation of Eq. ()2-16p . If we apply L{dt + 
drjdx + dp) to ■01 and separate out terms containing from the outcome, then 
from Eq. (|2-4p it must be identical to the coefficient of in 

[i(,dx + dp)Y{i;o-i^l). (2-18) 
Hence, explicitly written, Eq. ()2-16p reads 

r{dt + dr+ [lidx + dp)Y [{dx + dp)^ + k^o?]ARcPo 

= eV[i{dx + dp)Y{i>o-iPl) (2-19) 
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to 0(e). 

At this juncture, we note 

^PI-^q = B^ + 65|^|2 -B + 

|gife.=r^3^2^ + 3\A\'^A -A)+ e^^^ '^WA^ + q^^^-^A^ + c.c.} .(2-20) 

This, together with the identity: 

{d^ + dpff{p)^^ = MCk-kl)f[p), Ck = dl + 2ik-dp (2-21) 

for an arbitrary function of p, /(p), enable us to write down Eq. (|2-19p exphcitly. 
That is, 

drAR + ClAR = e{A-?,\A\''A) (NCOP) (2-22) 

for p = 0, and 

drAR-{Ck-kl)ClAR = e{Ck-kl){W''A + 'i\A\''A-A) (COP) (2-23) 
for p = 2. 

To calculate Eq. (j2-17p in the case p = 2, we may follow much the same procedure 
as that used above to derive Eqs. (I2^ and (f2^ from Eq. ([2461) . The 0(e) result 
is 

drBR - {dl + klfd^BR = ed^p{B^ + 6\A\^B - B). (2-24) 

Here we use A = ArZ, B = BrW on the RHS in Eqs. ([2^ - (l2^ . and then 
we may set p = and hence Z = W = 1 since the RHS should not depend on p 
explicitly. Therefore we finally obtain the proto-RG equations. Replacing the now 
dummy variables as Ar — )• A, Br B, p ^ x,t ^ t, they read 

dtA = {e-a^)A-3e\A\^A (NCOP), (2-25) 

and 

dtA = (□ - k^) [(D^ - e)A + 3e(|^|2 + B^)A] , (2-26) 

dtB = V^[{V^ + k^fB + e{B''^-B) + 6e\A\^B] (COP), (2-27) 

where □ is the rotationally covariant operator 

□ = V2 + 2ife-V. (2-28) 

The result ()2-25p is exactly of the same form as the amplitude equation for the 
SH model that Gunaratne et al^^ obtained first with a use of the multiple-scales 
method. Naturally, it also agrees with the proto-RG equation that Oono et al^ de- 
rived when they first formulated the proto-RG scheme. A set of proto-RG equations 
(f2^ and ([2^ is new. As a matter of fact, Eqs. ([2^ and ([2^ are identical 
to the results given in Ref. 112}) which are obtained with the multiple-scales method 
truncated at 0(e''/^); see, however, the caveat discussed in the next section. 
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§3. Discussion 

Notice that the appearance of the operator □ — A;q in Eq. (j2-26|) can be directly 
traced back to the use of the operator (dk + dpf' in the formula ()2-16p . In fact, each 
of the higher-order contributions from the RHS of this formula adds to the RHS of 
Eq. ([2^ the terms of the form - fcg) x function of ^, S, n = 2, 3, • • • . Thus 
Eq. (|2-16p in the case of COP yields the proto-RG equation of the form 

dtA = (□ - kl) [(rf - e)A + ^e{\A\^ + B^)A + h.o.t] , (3-1) 

where h.o.t. denotes the higher-order (in e) terms. Multiplying by and integrating 
the above equation ()3-ip . we find 



— / dxA{t,x)<Po{x) = (3-2) 



owing to the presence of the operator □ — /cq. The similar line of argument can be 
applied to Eq. (|2-17p . Thus we find the proto-RG equation for B takes the form 
dtB = V^[- • • ], which then guarantees that 

— / dxB{t,x) = 0. (3-3) 

Taken altogether, Eqs. ()3-2p and ()3-3p are in accord with the order-parameter conser- 
vation, J (1x1^(1, x) =constant, inherent in the COP equation. This in turn assures 
that the proto-RG equation respects the conservation law (and the rotational invari- 
ance as well) irrespective of at which order in e the perturbation expansion in the 
RHS of Eqs. OTUD and (I^TTj) is truncated. 

On the other hand, it is not the case with the multiple-scales formalism. The fact 
that in this formalism the conserving operator O — kQ is of mixed-order in e forces the 
very important conclusion that the finite-order result always breaks the conservation 
law EH To see this, let us recall the basic strategy of the so-called multiple-scales 
analysis; the reader should be warned that in fact it is not a multiple-scales (MS) 
analysis alone that is used but the combination of the MS analysis and reconstitution 
(hence called an MSR method). In the MSR analysis, first we carry out the MS 
analysis. Namely, we expand ■0 as a power series in the parameter S where e = 5^; 

tP = xPo + 6ipi + (5^2 + • • • , (3-4) 

and introduce slow space and time variables: 

X = 6x, T = 6^t. (3-5) 

The solution can be written as 

= (Aio(r, X) + 5A2o{T, X) + • • • ) + 

{ {Au (T, X) + 5A21 (T, X) + • • • ) e^*'-" + • • • + c.c} , (3-6) 

*-* Incidentally, the importance of realizing that the rotationally covariant operator □ is of mixed- 
order in e in the multiple-scales analysis for the SH model is emphasized in Ref. Ilip. 
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where the last ehipsis represents terms whose wavenumber differs from critical, being 
generated by the nonlinear interactions between the basic modes. At respective 
orders in 5, the Predholm alternative imposes the condition on the dynamics of 
Aii,AiQ,---. We now employ a kind of resummation procedure to combine the 
obtained conditions. To that end we introduce the reconstituted amplitude functions 
for the critical mode: 

A = Au + 6A21 + ■■■ , (3-7) 

as well as for the zero mode: B = Aiq + 5^20 + • • • • We can then deduce the 
amplitude equations for the dynamics of A and B. The operator ()2-28p is the spatial 
derivatives acting on the envelope function A and operates in the combination 

2i6k ■ Vx + -J^V^ 

where Vx operates on the slow variables X. When one generates the amplitude 
equation through the MSR method and if its perturbative expansion is carried out 
to all orders, then the amplitude equation will certainly guarantee the conservation 
law of the order-parameter field. However, this is not the case for finite truncation 
of the perturbation. For example, the correct term proportional to 

e\D-kl)A\A\' 

will appear on the RHS of Eq. (|3-ip if the amplitude equation is truncated at 0(5^), 
whereas at O(J^) we only obtain — e'^A;Q^|>l|'^ and with this term present the con- 
servation is lost. To state it differently, the symmetry-breaking term — e^feQ^I^I^ 
appearing at lower order in the perturbation expansion is symmetrized to restore 
the conservative symmetry only when the expansion is carried out to higher order. 
Since at any finite order there always appear such non-conserving terms, we cannot 
escape the afore-mentioned conclusion. Thus we see that, in the MSR analysis, in 
order to get the amplitude equation equivalent to our RG form (j2-26p one needs to 
"enforce" by hand the required conservation, which is obviously not a systematic 
step in the perturbation theory and is an uncontrollable approximation. In this con- 
nection one must always keep in mind that the MSR method is sometimes unreliable, 
and we refer the reader to Ref. IT3|) for the danger of the reconstitution. 

§4. Summary and Remarks 

We have presented a proper version of the proto-RG scheme to derive amplitude 
equations from NCOP and COP equations that produce spatially periodic patterns. 
The coarse-graining of the order parameter equations was performed by the proto-RG 
method. Our method is generically implemented the way in which the conservation 
law is preserved in the case of COP. There are two important aspects to be noted in 
our present formulation. 

The first is that 'integral constants' B and A in the solution (|2-7p depend on 
the renormalization point p. This lends a contrast to the conventional proto-RG 
formalism.-J The underlying idea is akin to the one imparted in the geometrical 
interpretation of the RG methocP^ in terms of the theory of envelopes. It was shown 
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there (we are paraphrasing here in terms of envelope curves) that the RG method 
is a theory manipulating the 'boundary conditions' at x = p. That is, suppose we 
have a perturbative solution which is valid only locally around x = p. Taking then 
p as a running parameter, we have a family of curves (i.e., perturbative solutions), 
each curve being a good approximation around x = p. Then it turns out that the 
envelope which contacts with each local curve at the point of tangency x = p gives 
a renormalized (global) solution. 

Secondly, Eqs. (|2-25p - (|2-27p are the 0(e) formal calculation result of the proto- 
RG formulation. The order of the operator □ — /cq or □ is not yet determined in this 
method. In order to complete the RG procedure, we must carry out the reduction 
further by introducing the multiple scales like (j3-5p . The resulting equation is called 
the RG equation.'^ Basically all that is required of such a reduction from proto-RG 
to RG equation is to find out the appropriate e-dependence of space-time scales, or 
what phenomenon you wish to see, which makes the RG equation consistent. In this 
reduction, there is always a risk of losing the relevant symmetry (and/or conservation 
law) for the reason expounded in §3. Your wish may or may not be granted. If you 
find the resulting RG equation is not consistent, so be it; it simply says your way 
of looking at the system is not consistent with the symmetry or conservation to 
the order you wish. For example, if you wish to require conservation to all orders 
(rigorous conservation) in the case of the proto-RG Eq. ()2-26p . with the choice of 
scales V ~ e^/^ ~ fc it is consistent with dt ~ e^/^. In this case, we have proto- 
RG eq. = RG eq., and consequently our RG equation preserves the symmetry and 
conservation law at the nontrivial lowest order. Finally, we remark that there is 
no generic reason why the RG equation should agree with an MSR result. It is 
particularly the case if the MSR result in hand is obtained with an inappropriate 
choice of the scales. A more detailed account and application of these ideas, however, 
is beyond the scope of the present paper, and it will be published in a separate paper. 
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